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The N-terminal region of influenza virus polymerase PB1 
adjacent to the PA binding site is involved in replication but 
not transcription of the viral genome 
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The influenza virus genome forms viral ribonucleoprotein (vRNP) complexes with 
nucleoprotein and viral RNA-dependent RNA polymerases (RdRp), PB1, PB2, and 
PA subunits. The vRNP complex catalyzes both genome replication and transcription 
reactions. PB1 contains the motifs highly conserved among RdRps and functions as 
a catalytic subunit of RdRp. The N-terminal region of PB1 between amino acid (a. a.) 
positions 1-83 contains both putative vRNA and cRNA promoter binding sites and a PA 
binding site. However, except for the PA binding site, the crystal structure and the function 
of the N-terminal region of PB1 are poorly understood. Here, we have examined the 
functional structure of the N-terminal region of PB1. The regions between a. a. positions 
1-50 are highly conserved between influenza A and B viruses, but amino acids at 
positions 16, 27, and 44 are different between two viruses. To elucidate the functional 
importance of these amino acids in replication and transcription of the viral genome, 
we generated viruses containing mutations at these positions by reverse genetics and 
examined replication and transcription activities of these mutants. We found that a. a. 
positions 27 and 44 are responsible for the viral replication activity but not transcription 
activity. 
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INTRODUCTION 

Influenza A and B viruses contain eight- segmented and negative - 
stranded RNAs (vRNA) as its genome. Each segment is 
encapsidated by nucleoprotein (NP) and associated with viral 
RNA-dependent RNA polymerases (RdRp) to form viral ribonu- 
cleoprotein (vRNP) complexes. The vRNP complex is a basic 
unit for both genome replication and transcription (Nagata et al, 
2008). 

The viral RdRp is a heterotrimer consisting of PB1, PB2, and 
PA subunits. Among them, PB1 functions as a catalytic subunit 
and assembly core of RdRp (Biswas and Nayak, 1996; Gonzalez 
et al, 1996; Toyoda et al, 1996; Zurcher et al, 1996; Ohtsu 
et al., 2002). The crystal structure of the interaction domains 
between N-terminal region of PB1 and C-terminal region of PA, 
and between C-terminal region of PB1 and N-terminal region 
of PB2 were resolved (He et al, 2008; Obayashi et al, 2008). 
PB1 contains the motifs highly conserved among RdRps, putative 
nucleotide-binding sites, and vRNA and cRNA promoter bind- 
ing sites (Asano and Ishihama, 1997; Li et al, 1998; Gonzalez and 
Ortin, 1999a,b; Kolpashchikov et al, 2004) (Figure 1A). 

The N-terminal region of PB1 (1-83 a.a.) contains both puta- 
tive vRNA and cRNA promoter binding sites. However, except 
for the PA binding site (1-15 a.a.), the function of this region 
was poorly understood. Alignment of amino acid sequences 
revealved that the a.a. positions 1-50 was highly conserved 
between influenza A and B viruses except for the amino acid posi- 
tions 16, 27, and 44. To identify the functional importance of 



these positions for the viral RNA synthesis, we determined the 
replicational and transcriptional activities using mutant viruses. 
Our results strongly suggest that the a.a. positions 27 and 44 are 
involved in replication process but not transcription process. 

MATERIALS AND METHODS 
BIOLOGICAL MATERIALS 

Monolayer cultures of 293T and MDCK cells were main- 
tained at 37°C with 5% CO2 in Dulbecco's Modified Eagle 
Medium (DMEM) and minimal essential medium (MEM) 
(Nissui), respectively, supplemented with 10% fetal bovine 
serum (Bovogen). Influenza virus strain A/WSN/33 (WSN) 
was prepared as previously described (Kawaguchi et al., 2005). 
Cycloheximide (CHX) was purchased from Sigma- Aldrich. 

GENERATION OF RECOMBINANT VIRUSES 

To construct plasmids from which human DNA- dependent RNA 
polymerase I (Pol I) transcribes mutated vRNAs, we amplified 
fragments containing mutated segment 2 by PCR using a plasmid 
containing wild type WSN segments in pHH21 vector (Neumann 
et al., 1999) as a template with sets of phosphorylated primers 
(see Table SI in the supplemental material). The amplified PCR 
products were self-ligated followed by sequencing. To generate 
recombinant viruses containing viral RNAs of WSN and mutated 
segment 2, reverse genetics system was used as described pre- 
viously (Neumann et al., 1999). After 48 h post transfection 
(hpt), aliquots of cell culture supernatants were used for virus 
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FIGURE 1 | Structure of PB1 subunit. (A) A diagrammatic representation of 
PB1. Black bar, cRNA promoter binding sites; gray bar, vRNA promoter 
binding sites; black and white vertical stripes, nucleotide binding sites; black 
and white horizontal stripes, PB2 binding site; black waved lines, PA binding 
site, and black dots, RNA dependent-RNA polymerase catalytic domain. 
(B) Alignment of amino acid sequences of putative RNA binding region 
common to vRNA and cRNA promoters (1-83 a. a.) among influenza A and 
influenza B viruses. PB1 sequences of A/WSN/33, A/WSN/1933 (H1N1); 



A/England/67, A/England/10/67 (H2N2); A/Hong Kong/68, A/Hong 
Kong/1/1968 (H3N2); A/Narita/98, A/parakeet/Narita/92A/98 (H9N2). These 
strains are picked up from different clades and periods of PB1 gene 
phylogenetic tree (Taubenberger et al., 2005); B/Aichi/88, B/Aichi/5/88; 
B/Alaska/92, B/Alaska/03/1992; B/Argentina/01, B/Argentina/1 32/2001 ; and 
B/Hong Kong/93, B/Hong Kong/02/1993. Sequences were aligned with 
CLUSTAL W. Asterisk (*), identical residues; colon (:), conserved residues, 
and dot (.), semi-conserved residues. 



amplification in MDCK cells. At 48 h post infection (hpi), the cul- 
ture fluid was collected, and the virus titer of these recombinant 
viruses was determined by plaque assays. 

RNA ANALYSIS BY qRT-PCR 

MDCK cells were infected with recombinant viruses at 
the multiplicity of infection (MOI) of 2.5. At 9 hpi, total 
RNA was isolated by the acid guanidine-phenol-chloroform 
method. To measure the accumulation levels of viral mRNA, 
cRNA, and vRNA, quantitative RT-PCR (qRT-PCR) was per- 
formed. Total RNAs were subjected to reverse transcription 
using ReverTraAce (Toyobo) with either (i) oligo (dT)2o> 
(ii) 5' - AGTAGAAACAAGGGTATTTTTCTTTA- 3' , or (hi) 
S'-GACGATGCAACGGCTGGTCTG-S' for synthesizing cDNA 
from segment 5 mRNA, cRNA, and vRNA, respectively 
(Kawaguchi and Nagata, 2007; Sugiyama et al., 2009). The 
synthesized single-stranded cDNAs were subjected to real-time 
quantitative PCR analysis (Thermal Cycler Dice real-time system 
TP800; TaKaRa) with SYBR Premix Ex Taq (TaKaRa) and a set 
of specific primers for segment 5 cDNA (see supplementary 
methods). The levels of these RNAs were normalized by the 



amount of cellular |3-actin mRNA measured using specific 
primers (see supplementary methods). These results are averages 
from three independent experiments with standard deviations. 
The level of significance was determined by Student's f-test 
(unpaired). 

RESULTS 

RNA SYNTHESIS OF INFLUENZA A MUTANT VIRUSES CONTAINING 
INFLUENZA B VIRUS-TYPE AMINO ACID SIGNATURES 

The N-terminal region of PB1 (1-83 a.a.) contains the PA bind- 
ing site and both putative vRNA and cRNA promoter bind- 
ing sites (Figure 1A). It is shown by alignment of amino acid 
sequences that the PB1 region between a.a. positions 1-50 are 
highly conserved between influenza A and B viruses, while the 
region between a.a. positions 51-83 differ between two viruses 
(Figure IB). In the highly conserved region, except for the PA 
binding site, a.a. at positions 16, 27, and 44 are different between 
these viruses. Furthermore, these a.a. positions in PB1 are con- 
served more than 99% of influenza A and B viruses deposited 
in the NCBI Influenza Virus Sequence Database (Table 1). To 
elucidate the functional importance of these a.a. for viral RNA 
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FIGURE 2 | Comparison of RNA synthesis of mutant influenza A 
viruses containing amino acids specific for influenza B viruses. The 

amino acids, N, D, and N at the positions 16, 27, and 44 were mutated to A, 
V, and I, respectively. MDCK cells were infected with mutant viruses at 
MOI of 2.5. At 9 hpi, the accumulation levels of viral mRNA (A), cRNA (B), 
and vRNA (C) were measured by qRT-PCR, and the amounts of these RNAs 
were normalized by that of cellular actin mRNA. These results are averages 
from three independent experiments with standard deviations. *p < 0.05; 
**p < 0.01. 



Table 1 | Conservation of amino acid position 16, 27, and 44 in PB1. 



Virus type a. a. position 
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The conservation of amino acid position 16, 27, and 44 in FBI was calcu- 
lated by using 7259 sequences of human and avian influenza A strains and 
1412 sequences of influenza B strains listed at the NCBI Influenza Sequence 
Database. 



synthesis, we generated influenza A viruses containing Ala at the 
a.a. position 16 (N16A), Val at the a.a. position 27 (D27V), and 
He at the a.a. position 44 (N44I) by reverse genetics. We exam- 
ined the RdRp activity by measuring the accumulation levels of 
viral mRNA, cRNA, and vRNA by qRT-PCR (Figure 2). The lev- 
els of all three type RNAs from D27V were increased compared 
with those from wild type and N16A virus, while those from N44I 
were significantly decreased. Based on the result that mutations 
at the positions 27 and 44 affect the synthesis activity of all viral 
RNAs equally, there could be two possibilities: (i) these muta- 
tions affect on the vRNA promoter recognition and followed by 
cRNA/mRNA synthesis, but do not affect on the cRNA promoter 
recognition and followed by vRNA synthesis, or (ii) these muta- 
tions affect independently the synthesis of each viral RNA, but 
total effects leads similar outputs in the synthesis of all viral RNAs. 

To elucidate whether these mutations affect genome repli- 
cation (cRNA and vRNA synthesis) and/or transcription (viral 
mRNA synthesis) activities, we measured the primary transcrip- 
tion activity using cycloheximide (CHX), a potent inhibitor of 
protein synthesis (Figure 3). It is shown that CHX suppresses 
viral protein synthesis and thereby leads to degradation of repli- 
cated virus genome RNA but not viral mRNA since newly vRNP 
formation was repressed (Vreede et al., 2004; Kawaguchi et al., 
2005). We utilized this method to measure the primary transcrip- 
tion activity that depends just only on incoming vRNP and is 



not affected by the replication process. In the presence of CHX, 
the levels of viral mRNA and vRNA were measured by qRT-PCR, 
and the transcription activity was represented as a ratio of viral 
mRNA/vRNA. This result shows that the transcription activity 
is not affected by these mutations, and thereby strongly suggests 
that these mutations affect the replication activity. 

AMINO ACID PROPERTIES AT a.a. POSITIONS 27 AND 44 FOR THE RNA 
SYNTHESIS ACTIVITY 

Aspartate at the position 27 is highly conserved among influenza 
A viruses, except for an H4N8 strain isolated from least sandpiper 
that contains asparagine (GenBank: ACI90 144.1). We generated 
D27E and D27N in addition to D27V (Figure 4). The RNA levels 
of mRNA, cRNA, and vRNA of D27N and D27V were increased 
significantly compared with those of wild-type and D27E. One 
of possible interpretations is that uncharged amino acids at a.a. 
position 27 may enhance the RNA synthesis. 

The mutation at the a.a. position 44 reduced the replica- 
tion activity (Figure 2). To clarify the importance of the a.a. at 
this position, we additionally generated N44D and N44Q viruses 
in addition to N44I and examined the RNA synthesis activity 
(Figure 5). The synthesis level of each viral RNA of N44I was 
decreased largely, while the amounts of mRNA and cRNA of 
N44D and N44Q were similar to those of wild type. In addition, 
the amount of vRNA of N44Q was more than that of wild type. 
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FIGURE 3 | Primary transcription activity of mutant influenza A viruses 
containing amino acids specific to influenza B viruses. MDCK cells 
were infected with mutant viruses at MOI of 2.5 and incubated in the 
presence of 1.0 |xg/ml of CHX. At 9 hpi, the accumulation levels of viral 
mRNA and vRNA were measured by qRT-PCR, and the amounts of these 
RNAs were normalized by that of cellular actin mRNA. The transcription 
activity is represented as a ratio of the amount of viral mRNA to that of 
vRNA. These results are averages from three independent experiments 
with standard deviations. 
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FIGURE 4 | RNA synthesis of viruses containing mutations at the 
amino acid position 27. Wild-type virus and mutant viruses containing 
amino acids D and E, N, and V, respectively, at the a. a. position 27 were 
infected into MDCK cells at MOI of 2.5. At 9 hpi, the accumulation levels of 
viral mRNA (A), cRNA (B), and vRNA (C) were measured by qRT-PCR, and 
the amounts of these RNAs were normalized by that of cellular actin 
mRNA. These results are averages from three independent experiments 
with standard deviations. *p < 0.05; **p < 0.01 . 



Thus, it is expected that the a.a. position 44 might be a water- 
soluble characteristic, and especially glutamine at this position 
stimulates the vRNA synthesis. 

DISCUSSION 

In this report, we have studied on three a.a. positions, i. e., 16, 
27, and 44, which are not conserved between influenza A and B 
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FIGURE 5 | RNA synthesis of viruses containing mutations at the 
amino acid position 44. Wild-type virus and mutant viruses containing 
amino acids N and I, D, and Q, respectively, at the a.a. position 44 were 
infected into MDCK cells at MOI of 2.5. At 9 hpi, the accumulation levels of 
viral mRNA (A), cRNA (B), and vRNA (C) were measured by qRT-PCR, and 
the amounts of these RNAs were normalized by that of cellular actin 
mRNA. These results are averages from three independent experiments 
with standard deviations. *p < 0.05; **p < 0.01 . 



viruses. The RNA synthesis activity of D27V was enhanced, while 
that of N44I was decreased (Figure 2). Based on these, we carried 
out further mutational analyses. The N44I showed the decreased 
level of RNA synthesis in three types of viral RNAs, while N44D 
did not affect the RNA synthesis (Figure 5). Interestingly, N44Q 
increased vRNA synthesis with little effect on viral mRNA and 
cRNA synthesis. It is possible that side chain group of Q may 
stimulate the cRNA promoter binding and increase the vRNA 
synthesis activity. 

D27V and D27N increased the RNA synthesis, while D27E 
mutation gave no effects (Figure 4). Furthermore, when the 
amounts of RNAs of D27V were analyzed at various MOI, those 
of D27V were increased (Figure S3). Although uncharged amino 
acid at this position enhances the RNA synthesis, molecular evo- 
lution has selected negatively charged amino acids. Therefore, it 
is assumed that charged amino acids at this position, even with 
low efficiency for the replication, are needed for PB1. Recently, 
mutational analyses showed that the sequences surrounding the 
PB1 AUG codon are multifunctional, and contain overlapping 
signals for translation initiation and for segment specific packag- 
ing (Wise et al., 201 1). We may consider a possibility that there is 
some regulatory coupling between replication and packaging and 
the a.a. position 27 has a role in this hypothetical mechanism. 

These a.a. positions are close to the PA binding site, and the 
PB1-RNA interaction could be affected by the presence of PA. 
We examined whether these mutations affect the assembly of 
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RdRp (Figure SI). The assembly of PB1 with PA and PB2 was 
not affected by these mutations. Moreover, these mutations did 
not affect the transcription activity, mRNA synthesis from vRNA 
(Figure 3 and Figure S2). Taken altogether, it is quite likely that 
amino acids at the positions 27 and 44 are involved in the repli- 
cation activity, possibly in cRNA promoter recognition with little 
effects on the transcription activity and the assembly of the RdRp 
complex. 

Recognition of the vRNA promoter depends on the 5 r -arm 
of the promoter and this binding improves the weak binding 
of RdRp to the 3 r -arm of the vRNA promoter (Tiley et al., 
1994; Gonzalez and Ortin, 1999a; Jung and Brownlee, 2006). 
Recognition of the cRNA promoter by RdRp has been shown by 
the in vitro binding of the PB1 subunit with the b'- and Smarms 
of the cRNA promoter (Gonzalez and Ortin, 1999b). Flexibility 
within the two uridines of the internal loop of the cRNA pro- 
moter required for protein binding in the cRNP complex (Park 
et al, 2003). Biochemical studies have shown that conformational 
changes in PB 1 of the influenza A virus RdRp lead to the inter- 
action with either vRNA or cRNA (Gonzalez and Ortin, 1999b). 
Based on previous reports and our findings, the positions 27 and 
44 may affect the PB1 structure, resulting in affecting PB1 binding 
activity to the 3 7 -arm of the cRNA promoter. Thus, we would pro- 
pose that these positions in PB1 are important for the replication 
activity by recognizing the cRNA promoter. 

ACKNOWLEDGMENTS 

This research was supported in part by a grant-in-aid from the 
Ministry of Education, Culture, Sports, Science, and Technology 
of Japan (to Kyosuke Nagata). 

SUPPLEMENTARY MATERIAL 

The Supplementary Material for this article can be found online 
at: http : //www. frontier sin . or g/j ournal/ 10.33 8 9/fmicb. 2 013. 
00398/abstract 

REFERENCES 

Asano, Y., and Ishihama, A. (1997). Identification of two nucleotide-binding 
domains on the PB1 subunit of influenza virus RNA polymerase. /. Biochem. 
122, 627-634. doi: 10.1093/oxfordjournals.jbchem.a021799 

Biswas, S. K., and Nayak, D. P. (1996). Influenza virus polymerase basic protein 
1 interacts with influenza virus polymerase basic protein 2 at multiple sites. 
/. Virol. 70, 6716-6722. 

Gonzalez, S., and Ortin, J. (1999a). Characterization of influenza virus PB1 pro- 
tein binding to viral RNA: two separate regions of the protein contribute to the 
interaction domain. /. Virol. 73, 631-637. 

Gonzalez, S., and Ortin, J. (1999b). Distinct regions of influenza virus PB1 poly- 
merase subunit recognize vRNA and cRNA templates. EMBO J. 18, 3767-3775. 
doi: 10. 1093/emboj/18. 13.3767 

Gonzalez, S., Zurcher, T., and Ortin, J. (1996). Identification of two separate 
domains in the influenza virus PB1 protein involved in the interac- 
tion with the PB2 and PA subunits: a model for the viral RNA poly- 
merase structure. Nucleic Acids Res. 24, 4456-4463. doi: 10.1093/nar/24.2 
2.4456 

He, X., Zhou, J., Bartlam, M., Zhang, R., Ma, J., Lou, Z., et al. (2008). Crystal struc- 
ture of the polymerase PA(C)-PB1(N) complex from an avian influenza H5N1 
virus. Nature 454, 1123-1126. doi: 10.1038/nature07120 

Jung, T. E., and Brownlee, G. G. (2006). A new promoter-binding site in the PB1 
subunit of the influenza A virus polymerase. /. Gen. Virol. 87, 679-688. doi: 
10.1099/vir.0.81453-0 

Kawaguchi, A., and Nagata, K. (2007). De novo replication of the influenza virus 
RNA genome is regulated by DNA replicative helicase, MCM. EMBO J. 26, 
4566-4575. doi: 10.1038/sj.emboj.7601881 



Kawaguchi, A., Naito, T., and Nagata, K. (2005). Involvement of influenza virus 
PA subunit in assembly of functional RNA polymerase complexes. /. Virol. 79, 
732-744. doi: 10.1128/JVI.79.2.732-744.2005 

Kolpashchikov, D. M., Honda, A., and Ishihama, A. (2004). Structure -function rela- 
tionship of the influenza virus RNA polymerase: primer-binding site on the PB 1 
subunit. Biochemistry 43, 5882-5887. doi: 10.1021/bi036139e 

Li, M. L., Ramirez, B. C, and Krug, R. M. (1998). RNA-dependent activation of 
primer RNA production by influenza virus polymerase: different regions of the 
same protein subunit constitute the two required RNA-binding sites. EMBO J. 
17, 5844-5852. doi: 10. 1093/emboj/17. 19.5844 

Nagata, K., Kawaguchi, A., and Naito, T. (2008). Host factors for replication and 
transcription of the influenza virus genome. Rev. Med. Virol. 18, 247-260. doi: 
10.1002/rmv.575 

Neumann, G., Watanabe, T., Ito, H., Watanabe, S., Goto, H., Gao, P., et al. (1999). 
Generation of influenza A viruses entirely from cloned cDNAs. Proc. Natl. Acad. 
Sci. U.S.A. 96, 9345-9350. doi: 10. 1073/pnas.96. 16.9345 

Obayashi, E., Yoshida, H., Kawai, E, Shibayama, N., Kawaguchi, A., Nagata, K., et al. 
(2008). The structural basis for an essential subunit interaction in influenza 
virus RNA polymerase. Nature 454, 1127-1131. doi: 10.1038/nature07225 

Ohtsu, Y., Honda, Y., Sakata, Y., Kato, H., and Toyoda, T. (2002). Fine mapping 
of the subunit binding sites of influenza virus RNA polymerase. Microbiol. 
Immunol. 46, 167-175. doi: 10.1111/j.l348-0421.2002.tb02682.x 

Park, C. J., Bae, S. H., Lee, M. K., Varani, G., and Choi, B. S. (2003). Solution 
structure of the influenza A virus cRNA promoter: implications for differential 
recognition of viral promoter structures by RNA-dependent RNA polymerase. 
Nucleic Acids Res. 31, 2824-2832. doi: 10.1093/nar/gkg387 

Sugiyama, K., Obayashi, E., Kawaguchi, A., Suzuki, Y., Tame, J. R., Nagata, 
K., et al. (2009). Structural insight into the essential PB1-PB2 subunit con- 
tact of the influenza virus RNA polymerase. EMBO J. 28, 1803-1811. doi: 
10.1038/emboj.2009.138 

Taubenberger, J. K., Reid, A. H., Lourens, R. M., Wang, R., Jin, G., and Fanning, 
T. G. (2005). Characterization of the 1918 influenza virus polymerase genes. 
Nature 437, 889-893. doi: 10.1038/nature04230 

Tiley, L. S., Hagen, M., Matthews, J. T., and Krystal, M. (1994). Sequence-specific 
binding of the influenza virus RNA polymerase to sequences located at the 5' 
ends of the viral RNAs. /. Virol. 68, 5108-51 16. 

Toyoda, T., Adyshev, D. M., Kobayashi, M., Iwata, A., and Ishihama, A. (1996). 
Molecular assembly of the influenza virus RNA polymerase: determination 
of the subunit- subunit contact sites. /. Gen. Virol. 77 (pt 9), 2149-2157. doi: 
10. 1099/0022- 13 17-77-9-2 149 

Vreede, F. T., Jung, T. E., and Brownlee, G. G. (2004). Model suggesting that replica- 
tion of influenza virus is regulated by stabilization of replicative intermediates. 
/. Virol. 78, 9568-9572. doi: 10.1128/JVI.78.17.9568-9572.2004 

Wise, H. M., Barbezange, C, Jagger, B. W, Dalton, R. M., Gog, J. R., Curran, 
M. D., et al. (2011). Overlapping signals for translational regulation and pack- 
aging of influenza A virus segment 2. Nucleic Acids Res. 39, 7775-7790. doi: 
10.1093/nar/gkr487 

Zurcher, T., de La Luna, S., Sanz-Ezquerro, J. J., Nieto, A., and Ortin, J. (1996). 
Mutational analysis of the influenza virus A/Victoria/3/75 PA protein: studies of 
interaction with PB1 protein and identification of a dominant negative mutant. 
/. Gen. Virol. 77 (pt 8), 1745-1749. doi: 10.1099/0022-1317-77-8-1745 

Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that could be 
construed as a potential conflict of interest. 

Received: 17 October 2013; accepted: 04 December 2013; published online: 18 
December 2013. 

Citation: Binh NT, Wakai C, Kawaguchi A and Nagata K (2013) The N-terminal 
region of influenza virus polymerase PB1 adjacent to the PA binding site is involved 
in replication but not transcription of the viral genome. Front. Microbiol. 4:398. doi: 
1 0. 3389/ fmicb. 201 3. 00398 

This article was submitted to Virology, a section of the journal Frontiers in 
Microbiology. 

Copyright © 2013 Binh, Wakai, Kawaguchi and Nagata. This is an open-access arti- 
cle distributed under the terms of the Creative Commons Attribution License (CC BY). 
The use, distribution or reproduction in other forums is permitted, provided the 
original author(s) or licensor are credited and that the original publication in this 
journal is cited, in accordance with accepted academic practice. No use, distribution or 
reproduction is permitted which does not comply with these terms. 



www.frontiersin.org 



December 2013 | Volume 4 | Article 398 | 5 



